Geometry- and diffraction-independent ionization probabilities in intense laser fields: 
probing atomic ionization mechanisms with effective intensity matching 
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We report a novel experimental technique for the comparison of ionization processes in ultrafast 
laser pulses irrespective of pulse ellipticity. Multiple ionization of xenon by 50 fs 790 nm, linearly 
and circularly polarized laser pulses is observed over the intensity range 10 TW/cm 2 to 10 PW/cm 2 
using Effective Intensity Matching (EIM), which is coupled with Intensity Selective Scanning (ISS) 
to recover the geometry-independent probability of ionization. Such measurements, made possible 
by quantifying diffraction effects in the laser focus, are compared directly to theoretical predictions of 
multiphoton, tunnel and field ionization, and a remarkable agreement demonstrated. EIM-ISS allows 
the straightforward quantification of the probability of recollision ionization in a linearly polarized 
laser pulse. Furthermore, probability of ionization is discussed in terms of the Keldysh adiabaticity 
parameter 7, and the influence of the precursor ionic states present in recollision ionization is 
observed for the first time. 

PACS numbers: 32.80.Rm, 42.50.Hz, 42.65.Sf 



I. INTRODUCTION 



Over the last decade and a half, the advent of modern 
laser technology has resulted in a fruition of significant 
experimental studies on the interaction of intense ultra- 
fast laser pulses with dilute matter. The success of these 
works has naturally precipitated a reaction from the the- 
oretical community, which has lead to a highly productive 
dialogue. Particular highlights are the observation p| of 
high harmonic generation (HHG), the use of HHG to gen- 
erate attosecond-timescale bursts of XUV radiation 0, 
nonsequential ionization of atoms Q and molecules |j 
through electron rescattering and the enhanced ioniza- 
tion of molecules leading to Coulomb explosion [|| . Ex- 
perimental advances by the authors and co-workers have 
allowed the investigation of the interaction of ultrafast 
laser pulses with ionic targets, where the ionization of 
ground and metastable Ar + C + Q and Xe + Q has 
been observed for the first time. Central to all of these 
studies is the observation of process which are depen- 
dent on focused laser intensity, traditionally facilitated 
by variation of the pulse energy. While being straight- 
forward to carry out experimentally, this technique suf- 
fers the drawback that, as the pulse energy is varied, 
the size of the laser focus is altered. Such geometry- 
dependent measurements are of course still of great im- 
portance, and have allowed major advances in the under- 
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standing of laser-matter interactions: the drawback is in 
making comparisons with theory. A different approach 
has been facilitated by the ongoing development of laser 
amplification techniques 9] , permitting macroscopic con- 
trol of laser focal conditions through intensity selective 
scanning (ISS) 0, [n|. Here the laser focal geometry 
is constant, and a spatially selective detector is used to 
image different regions of the focus. Such measurements 
are of particular interest to the atomic physics commu- 
nity, with regards, for example, to the response of noble 
gas atoms to focused femtosecond laser pulses. 

While a number of experimental studies have com- 
pared atomic ionization using circular and linear polar- 
ized laser pulses, for recent examples see 0,^2, to the 
authors' knowledge, all have employed the traditional 
geometry-dependent technique. For a particular pulse 
energy, there is a change in electric field amplitude when 
the polarization is switched from linear to circular, mak- 
ing a comparison problematic. We propose a new method 
to circumvent this problem, as a development of ISS. A 
constant pulse energy ratio between linear and circular 
polarizations is defined such that the relative ionization 
yield is constant as the laser focus is translated past the 
spatially selective detector, irrespective of laser polar- 
ization. This technique of effective intensity matching 
(EIM) coupled with ISS is employed to examine ultra- 
fast strong field ionization of xenon: EIM-ISS results are 
presented. 

Before quantitative analysis of the EIM-ISS results is 
possible, it is necessary to remove the geometry depen- 
dence inherent in the ionization yields. A method to per- 
form such a deconvolution has been proposed in • An 



analogous technique has been presented by the authors 
|13|. where a numerical inversion is employed to remove 
the geometry dependence from the ISS results. However 
both techniques only allow the removal of Gaussian fo- 
cusing. Modern laser systems typically generate pulses 
with a Gaussian profile in the far-field. However, optical 
transport systems impose geometrical constraints such 
that it is practically impossible to propagate this per- 
fect profile into the experimental chamber. Such beam 
propagation has been addressed recently by Lii and co- 
workers 0| through the application of an aperture func- 
tion approximation |l5j to a multi-apertured optical sys- 
tem. Unfortunately, this approximation, even with little 
or no beam truncation, will introduce oscillatory distor- 
tions to the laser intensity. Zhang and co-workers |Tfij 
have presented an elegant analysis of the diffraction of 
a focused beam, with particular application to tradi- 
tional geometry-dependent intensity variation measure- 
ments. In the present work, beam diffraction is quan- 
tified before focusing takes place: this treatment better 
representing typical experimental systems. 

The spatial insensitivity of most instruments employed 
in intensity variation measurements 0, 0, 0, 0, [2(], 



l21l l22fl average over such diffraction effects, as all of the 
confocal volume lies within the volume that the instru- 
ment images. In the case of EIM-ISS, the very spatial 
sensitivity that makes the technique so powerful also 
makes it susceptible to optical distortions. Given that 
the EIM-ISS technique naturally discriminates between 
low-intensity large volume and high-intensity small vol- 
ume processes, then before the laser - atom system under 
investigation can be fully understood, it is vital to quan- 
tify the spatial distribution of laser intensity. Herein, we 
present the results of an analytical treatment of the fo- 
cusing of a truncated Gaussian laser pulse through an 
arbitrary ABCD (where A to D are elements of the Sys- 
tem Matrix) optical system, the derivation of which is 
given in Appendix A. 

By matching the optical conditions in our ABCD 
model to those in our experimental focus, we then re- 
move the geometry dependence from the EIM-ISS data, 
revealing for the first time geometry-free and diffraction- 
free probabilities of ionization for both linear and circular 
polarizations, and compare them directly to theoretical 
predictions. 



II. APPLICATION TO LASER - DILUTE 
MATTER INTERACTIONS 



In the following section, the optical system employed is 
defined, then the xenon EIM-ISS data is presented. The 
apertured solution for an arbitrary ABCD optical sys- 
tem is then employed to remove the geometrical influence 
of the volume of the laser focus, resulting in geometry- 
independent atomic ionization probabilities. Ionization 
mechanisms arc then discussed in detail. Finally, the 
variation of geometry-independent ionization probability 




FIG. 1: Illustration of the cylindrically symmetric optical sys- 
tem employed in recent laser - dilute matter interactions. The 
aperture, radius a, lies in the input plane at z = 0, with radial 
co-ordinate n . The resulting focus lies in the output plane at 
2 = zi + 22, with radial co-ordinate r%. The focal length of 
the lens is /, and Zf = zi - f is position with respect to the 
focus. 



with the Keldysh parameter |23j is presented, allowing 
ionization mechanisms to be discussed in terms of the rel- 
ative frequencies of the tunneling electron and the laser 
field. 



A. Optical system 

In recent experimental studies published by the au- 
thors and co-workers [ESS Ei3 the optical system is as 
illustrated in Fig. ^ with a corresponding system matrix 



A B 
C D 



1 - (z 2 /f) Zl + (1- Zl /f) Z2 



-l/f 



i - (V/) 



(1) 



where z\ and z 2 correspond to the distance between the 
aperture and lens, and lens and focus, and / is the fo- 
cal length of the lens. We also define Zf — z 2 - / as the 
position parallel to the beam propagation direction with 
respect to the focus. The system matrix is derived using 
straightforward matrix optics, and is the transformation 
from the input plane (containing the aperture) to the 
output plane, through a translation, a refraction and a 
second translation to the focus. In the present work, the 
following values apply: z\ = 300 mm, z 2 — 250 mm, / 
= 250 mm and the aperture radius a = 11 mm. Given 
the wavelength, A = 790 nm, the only quantity which is 
unknown is the beam radius uj g before the aperture at 
which diffraction occurs. The fixed aperture, located at 
2=0, defines the reference plane for the System Matrix. 
As presented in Appendix A, we derive a solution for 
the laser intensity distribution in the vicinity of the fo- 
cus, which accounts for diffraction of the incoming laser 
beam at an aperture of finite diameter. Such a solution is 
particularly applicable to ultrafast (Ti: sapphire based) 
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physics, where the intensity distribution in the laser fo- 
cus is rarely Gaussian. Fig. |2fa)-(j) shows the result of 
simulating the focus of the current optical system for the 
fixed radius aperture and a range of beam radii. Initially, 
the beam radius lo 9 — 5 mm, and the focus created is al- 
most identical to the unapertured form described by Eq. 
(A6) in Appendix A. However, as the beam radius is in- 
creased, the focus is disturbed by the diffraction of the 
incoming beam at the aperture. As the beam radius be- 
comes comparable to the aperture radius, the focus shows 
pronounced lines of maximum and minimum intensity in 
both Zf and r%. The isointensity contours in Fig. Efa)- 
(j), separated by an order of magnitude, illustrate that as 
the beam radius is increased, the focal spot size increases 
along T2 contrary to the behaviour expected in the case 
of an unapertured Gaussian beam. 



B. Effective intensity matching (EIM) 
as applied to ISS 



The experimental apparatus employed in the present 
work has been described in detail ja 0, H, 0, EL so 
only a brief outline is given here. The 20 mJ, 790 nm 
50 fs output of the ASTRA Ti:sapphirc laser is trans- 
mission focused on to the target gas in the source region 
of the time-of-flight mass spectrometer. Given the high 
power of the laser used in our recent work, intensities in 
excess of 10 17 Wcm~ 2 are routinely generated in a long 
focus. Indeed, the intensity range required for ultrafast 
ionization rate studies, typically > 10 13 Wcm~ 2 is pro- 
duced « 10 mm from the centre of the focus. The 250 
/jm entrance aperture located in the time-of-flight mass 
spectrometer places a tight spatial limitation on the vol- 
ume of the focused laser beam which produces ion signal 
from the target gas, and by translating the focusing lens 
parallel to the z-axis, the spectrometer is exposed to dif- 
ferent regions of the focus. Typically, such ISS or z-scan 
experiments are performed by averaging the ion signal 
measured at each lens position, Zf. 

To explore the ionization mechanism in xenon we have 
developed a novel experimental technique, namely effec- 
tive intensity matching, combined with intensity selective 
scanning (EIM-ISS). The essence of EIM-ISS is to define 
a constant ratio, J?eim between the laser intensities of 
the linearly (Jii n ) and circularly (J C irc) polarized laser 
beams such that the spatial distribution of the ions de- 
tected for each polarization are the same for all z/ values 
assuming that nonsequential ionization processes are neg- 
ligible. Before recording the xenon data presented, the 
ionization of neon was observed with circularly and lin- 
early polarized radiation to define the ratio Reim- Neon 
is a good test gas as, of all the noble gases, it is least 
susceptible to nonsequential (recollision) ionization [25) . 
Our studies, not presented here reveal that J?eim = 
Jiin / Jcirc = 0.65 ± 0.02 gives a remarkable match in 
the Ne + ionization signal over a large (« 10 mm) range 
of Zf, equivalent to an intensity range 10 13 Wcnr 2 to 
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FIG. 2: (a) to (j) Simulated intensity distributions in loga- 
rithmic grayscale with isointensity contours, separated by an 
order of magnitude, for the optical system in Fig. 1. The 
aperture radius a = 11 mm is kept constant, while the beam 
radius uo g is varied from 5 to 20 mm. Diffraction is clearly 
apparent even when the aperture is considerably larger than 
the beam radius. 



10 16 Wcm -2 . A similar approach was reported in [2g] 
to explain threshold ionization intensities observed with 
linearly and circularly polarized light in a long-pulse (1 
ps) traditional intensity variation measurement. How- 
ever, the present study and that of -jj are of a higher 
precision, being the culmination of a systematic inves- 
tigation, and are unique in the area of ultrafast intense 
field interactions. 

The exact value of J?eim is determined by the influ- 
ence of the different laser polarizations and electric field 
amplitudes on the ionization process. In the linearly po- 
larized case, the sinusoidal laser electric field is modu- 
lated by the pulse envelope: the field amplitude oscillates 
under a typically sech 2 temporal profile. With circular 
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Position w.r.t. focus, z f (mm) 

FIG. 3: Raw EIM-ISS data for the ionization of xenon to 
Xe n+ (n = 1, 2, 3). The integrated xenon ion yield is recorded 
as function of focusing optic position with respect to the axis 
of the spectrometer, in Fig. 1 equivilant to Zf. Linear (solid 
symbols) and circular (open symbols) polarized laser pulses 
are employed, where EIM is used to define the distribution 
of ion signal with varying zj . The presence of nonsequential 
(recollision) ionization is apparent in the Xe 2+ and Xe 3+ sig- 
nal as signal enhancement at low intensity (large Zf) in the 
case of linear polarization. Ripples in the ion yield indicate 
the presence of diffraction. 



polarization the electric field is continually present, and 
the electric field direction rotates through 2n during the 
laser period, and the field amplitude typically takes a 
sech 2 profile. In the circularly polarized field, the projec- 
tion of the angular momentum imparted to the electron 
is conserved along the direction of beam propagation. 
However, in the linearly polarized field, momentum is 
imparted to the electron in the direction of electric field 
oscillation, i.e. perpendicular to the propagation axis. 
This has an important influence on tunnel ionization, as 
recently noted by Tulenko and Zon |27j in which the rate 
of tunnelling depends on the magnetic quantum number 
of the tunnelling electron: in |2g it was concluded that 
the ionization rate should be identical, irrespective of po- 
larization type. 

With -Reim defined, the spectrometer source region was 
filled with xenon to a pressure of 1 x mbar: a low 

pressure is used to avoid the influence of space-charge 



effects. All ions generated were averaged over 500 laser 
shots for each position along Zf. The observed z-scan is 
presented in Fig. [3J showing Xe™ + (n — 1, 2, 3) product 
ion yields, with the Xe 2+ and Xe + yields displaced verti- 
cally by 0.5 and 1.0 respectively for clarity. By limiting 
the field of view of the spectrometer, the ion signal from 
a particular charge state at each z-position is an integral 
over all T2 and a narrow range of Zf. This selectivity is 
apparent in the data in Fig.|3| Starting at the focus, each 
ion yield presented is observed to increase from near zero 
to a maximum, the position of the maximum moving to 
lower Zf with increasing charge state, n. Each successive 
charge state requires an increase in intensity to maximize 
its ionization yield, a common result of traditional inten- 
sity measurements El E3 El El El EH El Subtle rip- 
ples present in the data presented in Fig.|31are the result 
of diffraction in the ultrafast laser focus. Importantly, 
linear and circular time-of-flight spectra are recorded at 
each Zf position as the focusing optic is translated, rather 
than recording a linear scan followed by a circular scan. 
This ensures that conditions within both the spectrome- 
ter and laser focus undergo minimal variation. 

Xenon has a far greater propensity for ionization 
through recollision in a linear laser field than neon, giv- 
ing rise to the major differences between linearly and 
circularly polarized light observed in the EIM-ISS data 
in Fig. |3 The data recorded with circular polarization is 
recollision- free, as the probability of the first-ionized elec- 
tron returning to the ion in a circularly polarized field is 
negligible |29j. The 'signature' of recollision is clearly ap- 
parent in the Xe 2+ and Xe 3+ data as an enhancement of 
the integrated ion yield with linear polarization at po- 
sitions of z at larger distances from the focus than the 
maxima (for example, in Xe 2+ 4.5 mm < Zf < 8 mm). 
This enhancement is due to recollision ionization in the 
linear field being more efficient than sequential multiple 
ionization by the laser field at low intensities (i.e. large 
Zf with respect to the centre of the focus). 

The enhancement of ion yield by recollision ionization 
in Xe 2+ and Xe 3+ is directly responsible for the signifi- 
cant differences in the Xe + yield curves. This is a con- 
sequence of the conservation of confocal volume, necessi- 
tated by the following important points. The core result 
of EIM-ISS is that the effective laser intensity is identi- 
cal for linear and circular polarizations. Given that the 
number density is constant, the ionization yield from se- 
quential ionization of the atoms in the confocal volume is 
made to match for both polarizations over all Zf. How- 
ever, the volume of the laser focus must be conserved, 
therefore, we observe for the first time the true interplay 
between the volumes generating a charge state Xe™ + (n 
1 . 2, 3...), which totally contains and often overlaps with 
the volume generating all higher charge states. Therefore 
the suppression of integrated ion yield in the case of Xe + 
in the range 3 mm < Zf < 8 mm is due to the depopula- 
tion of the Xe + volume in the laser focus (when the laser 
is linearly polarized) through the observed mechanism of 
recollision ionization to higher charge states. An analo- 
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gous depletion is observed around the maximum of the 
Xe 2+ yield, caused by depopulation of the Xe 2+ volume 
by recollision ionization to Xe 3+ in linearly polarized ra- 
diation. This sequence is repeated for the higher charge 
states. By comparing the sum of ion yields over all charge 
states at each z for the two polarizations, a near exact 
match is observed, confirming conservation of focal vol- 
ume. This then is the major benefit of EIM-ISS, allowing 
direct comparisons to be made between atomic ionization 
mechanisms for the first time. 



C. Deconvolution of EIM-ISS results: 
volume-independent ionization probabilities 

To further develop our understanding of the interplay 
of ionization mechanisms as observed using EIM-ISS us- 
ing both linear and circular polarizations, we turn now 
to removing the geometry dependence of the ionization 
signal from the EIM-ISS data in Fig. [3] The aim of such 
analysis is the recovery of geometry-independent ioniza- 
tion probabilities. 

Van Woerkom and co-workers have previously estab- 
lished a technique for the removal of the dependence of 
ionization probability on focal volume lOJllj. where the 
ionization probability il(/), may be calculated from the 
Zf dependent ion signal S(zf ) via the on-axis intensity 
distribution I ax {zf) according to: 



ft(Iax(zf)) OC 



dl ax (zf)/dzf ) dzf 



[I ax (z f )S(z f )] (2) 



In the present work, we have refined this treatment, 
firstly to allow for the non-Gaussian nature of the fo- 
cussed beam and secondly removing the dependence of 
ionization probability on the intensity distribution within 
the laser focus, despite the presence of diffraction. The 
latter is achieved by applying Eq. (A9) in Appendix A 
to Eq. (2) where I ax = \U(r 2 =0, z f ) 2 \. 

Although the 1/e 2 beam radius is difficult to define 
accurately, it is of a comparable size to the final aper- 
ture radius. Furthermore, the ionization signal of the 
first three charge states is expected to saturate within 
the range of z/ (i.e. intensity) in the current experiment 
|2l| , thus the gradient of the ionization probability should 
tend to zero at high intensity for all change states, n. At 
intensities above saturation, it is reasonable to expect 
ionization probability to decrease, as ionization to higher 
charge states becomes increasingly more likely. However, 
if we define some saturation intensity /sat, which oc- 
curs at a corresponding on-axis distance 2 sat from the 
focus, at Zf < zsat (he. on-axis I > /sat), I sat is 
still present off-axis thus occupying a larger proportion 
of the volume than the on-axis contribution. The result 
of applying Eq. (2) is only valid at intensities < I sat- at 
intensities > /sat, we define the probability of ionization 
as unity. This volume- and diffraction-free intensity de- 
pendent quantity is referred to as the partial probability 




Intensity (Wcm" ) 

FIG. 4: EIM geometry-independent partial probability of ion- 
ization (PPI) for Xe' 1+ (n = 1, 2, 3) as function of intensity, 
measured for both linear (solid symbols) and circular (open 
symbols) polarizations. Increasing charge states are normal- 
ized to descending orders of magnitude for clarity. Three 
types of sequential ionization are identified: multiphoton (low 
intensity, constant gradient, long dashed line), tunnelling (in- 
termediate intensity, decreasing gradient, solid line) and field 
ionization (high intensity, vertical arrow). Recollision ioniza- 
tion (short dashed line) is also observed with linear polariza- 
tion. See text for discussion. 



of ionization (PPI). If the conserved probability of ion- 
ization (CPI) is required, this may be readily calculated 
according to CPI (n) = PPI (n) - PPI (n + 1), where n 
is charge state, as before. This technique will be applied 
in future publications. 

The geometry- and diffraction-independent PPIs are 
now recovered from the EIM-ISS data in Fig. 03 and are 
presented in Fig. QJ where the saturated (unity) ioniza- 
tion probabilities of the increasing charge states are nor- 
malized to descending orders of magnitude to aid visual 
presentation. The beam radius u> g is estimated from a 
measurement of the unfocused beam profile, and small 
adjustments (of the order of 0.1 mm) made to uj g until 
the PPI curves recovered do not exhibit rapid changes of 
gradient with intensity. The beam radius oj g — 10.25 mm 
used to recover these ionization probabilities is found to 
be not only consistent over all charge states observed in 
xenon (n — 1, 2, 3), but also for a wide range of atomic 



6 



(the other Noble gases) and molecular (H2, D2, N2, CO2) 
targets covered in our studies. 

To quantify the results presented in Fig. 4, the 
intensity-dependent probability of ionization to Xe™ + (ra 
= 1, 2, 3) by four ionization mechanisms is now discussed. 
Starting with the lowest intensity, the following mecha- 
nisms are apparent in Fig. 

(i) Multiphoton ionization (MPI) is predicted by 
lowest-order perturbation theory (LOPT) to vary accord- 
ing to I N 30], where N is the number of 790 nm (1.56 
eV) photons absorbed, thus on a log-log plot of proba- 
bility vs. intensity a constant gradient of N should exist. 
Increasing charge state ra requires an increasing number 
of photons. For example, Xe 2+ production requires an 
energy of 21.21 eV, therefore the absorption of at least 
fourteen 790 nm photons is required. The expected gra- 
dients for MPI of Xe"+ are 7.8 (ra = 1), 13.6 (n = 2) 
and 20.6 (ra = 3). For n = 1 and 2, the expected gradi- 
ents are shown in Fig. 4 as long dashed lines. For ra = 
3, there is no MPI observed within the sensitivity of the 
experiment, determined by the low target gas density. 

(ii) Tunnel ionization (TI), as described initially by 
Keldysh [23] and later refined by Popov and co-workers 
[28|. is the result of a bound electron tunnelling out of 
the atom through the laser-modified Coulomb potential. 
With increasing intensity, TI becomes more efficient, as 
the laser + Coulomb barrier width decreases. Repre- 
sented in Fig. 4 as solid lines, Keldysh theory predicts 
a tunnel rate, which we have converted to probability. 
Note that the amount of TI observed is not always well- 
predicted by tunnelling theory, as observed recently by 
Yamakawa and co-workers |31| who introduced scaling 
factors in order to match theory and experiment. In the 
present work, when the TI theory is shifted to improve 
the fit with the PPI, the factors used are made clear. 

(iii) At the highest intensities, classical field ionization 
(FI) dominates, as the electric field of the laser pulse is 
sufficient to rapidly suppress the Coulomb potential, al- 
lowing direct liberation of the electron. The intensity at 
which FI is predicted to occur is indicated by the verti- 
cal arrows in Fig. 4, calculated using the over-the-barrier 
model. The intensity at which FI dominates should co- 
incide with the saturation of TI. 

Ionization mechanisms (i) to (iii) are sequential in na- 
ture, and due to the EIM technique should occur at the 
same effective laser intensity for linear and circular po- 
larizations (evident from the overlap of the Xe + PPIs 
in Fig. 4). As the ionization potential increases with 
charge state, we expect a systematic increase in the in- 
tensity required to produce a certain ionization mecha- 
nism. Furthermore, as discussed with reference to Fig. 3, 
there is an enhancement of signal in a linearly polarized 
laser pulse at low intensities, thus we must also consider 
nonsequential ionization processes: 

(iv) Recollision ionization (RI) [29|. also referred to 
as nonsequential double or multiple ionization (NSDI or 
NSMI) , predominantly occurs due to the linear laser field 
driving a correlated liberated electron back to the par- 



ent ion, thus initiating secondary ionization 32]. This 
process has been the subject of a number of COLTRIMS 
studies [13. In Fig. 4, the enhancement of Xe 2+ and 
Xe 3+ by RI is clearly seen in the difference of PPI curves 
for linear and circularly polarized radiation. In Fig. 4, 
the presence of RI is indicated by the short dashed lines, 
generated by transferring the theoretical prediction of se- 
quential ionization of previous charge state(s) onto the 
linearly polarized data. For Xe 2+ RI, the combined MPI 
and TI to Xe+ PPI is fitted, whereas for Xe 3 + RI, the 
MPI and TI to Xe+ and the TI to Xe 2+ is fitted as there 
is more than one recollision mechanism possible resulting 
in triple ionization. 

By comparing the quality of fit of MPI, TI, FI and RI 
to the data in Fig. 4, we can directly determine how suc- 
cessfully the geometry- and diffraction-independent PPIs 
are predicted by theory. As the PPI measured with circu- 
lar polarization is not influenced by RI, the circular data 
is therefore a more direct comparison to the predicted 
MPI and TI response. 

For ionization to Xe + , MPI with the absorption of 
eight photons occurs at the lowest intensities (< 2.5 x 
10 13 Wcm" 2 ) present in Fig. 4, indicating that LOPT 
is applicable. As the laser intensity increases, the PPI 
response tends away from MPI: the TI prediction fits 
the data excellently (i.e. directly predicted by Keldysh 
theory), even at intensities greater than that where FI 
is predicted to be the dominant mechanism. The Xe + 
data clearly illustrates the intensity ranges over which 
the different ionization mechanisms apply. 

When the laser field generates Xe 2+ , the absorption of 
fourteen photons is required for MPI to proceed. As is 
apparent from the long dashed line, only the very low- 
est intensities are predicted by MPI. TI theory is even 
more successful in this, accurately predicting a three or- 
der of magnitude increase in the PPI to better than 10 
%. However, it is necessary to translate the data by a 
factor of 0.91 in intensity to achieve this fit (cf 0.85 from 
|3l| . albeit during a shorter duration laser pulse, and 
geometry-dependent ion yield was measured). As with ra 
= 1, FI is only reasonably accurate. In the case of the 
linearly polarized laser pulse, at low intensity, there is 
considerable RI present, as indicated by the short dashed 
line. The RI contribution follows clearly the shape of the 
theoretical prediction of sequential MPI and TI produc- 
tion (short dashed line) visually fitted to the Xe 2+ data. 
This is consistent with Xe + being the source atom for 
RI to Xe 2+ . Such a method has been used by a number 
of groups (for example, see |2l|) fitting RI in geometry- 
dependent ion yield measurements, and it appears to be 
successful here. 

LOPT does not apply at all in the case of Xe 3+ as the 
expected gradient of 20.6 is far too high to be supported 
by the data. However, TI theory is able to accurately pre- 
dict the PPI from low intensity right up to saturation, 
with the data translated by a factor of 0.85 in inten- 
sity, (cf 0.70 from [3l|), here well-defined by FI (vertical 
arrow). There are now three possible nonsequential RI 
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routes possible: (a) — > 1 3, (b) => 2 — > 3 and 
(c) 0^3, where — > indicates sequential ionization, =>■ 
indicates double ionization through RI, and ^ indicates 
triple ionization through RI. While there is strong evi- 
dence for mechanisms (a) and (b) as is clear from the 
two short dashed lines on Fig. 4, there is little evidence 
for the presence of mechanism (c). If it does occur, it is 
with a greatly suppressed PPI as compared to (a) and 
(b). The shape of the RI PPI for n — 3 is well described 
by a combination of the theoretical prediction of sequen- 
tial ionization to Xe + and Xe 2+ as denoted by the short 
dashed lines (a) and (b) in Fig. 4 respectively. 

The ease with which ionization mechanism may be de- 
termined from Fig. 0] is illustrated by considering what 
happens to a group of atoms in a circularly polarized laser 
field at an intensity of 2 x 10 14 Wcm -2 . Such an inten- 
sity will at least triply ionize the atom: all atoms will 
be ionized to Xe + by FI, between 40% and 50% of these 
ions will be further ionized to Xe 2+ by TI, and of these 
ions, 1% will undergo MPI to Xe 3+ . The situation with 
linear polarization is complicated by NSMI, however the 
contribution is readily calculable from Fig. 0] Through a 
combination of multiphoton, tunnel and field ionization 
theory, the response of xenon to the ultrafast laser field 
may be accurately quantified. 
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FIG. 5: PPI as a function of Keldysh parameter for Xe™ + (n 
= 1, 2, 3) for (a) circular and (b) linear polarization. The 
presence of multiphoton ionization (MPI) dominates at high 
7, whereas tunnel ionization (TI) is more prominent at 7 < 2. 
A five point running average (solid line) through the points 
in (a) is duplicated onto (b) (dashed line) to illustrate the 
sequential ionization PPI as a function of 7. Interestingly, 
the contribution from recollision ionization in the linear pulse 
tends to MPI/TI response of Xe + . 



D. Ionization mechanism 
and the Keldysh parameter 

The pioneering theoretical work of Keldysh |2^| estab- 
lished a rule for distinguishing between ionization mech- 
anisms in a laser field using the fact that the ionization is 
governed by the relative frequency of the laser field and 
the tunnel frequency of the electron. The ratio of these 
frequencies is defined as the adiabaticity, or Keldysh pa- 
rameter, 7 [23ll28j|. and allows the ionization mechanism 
to be broadly determined: 



7 



^las 



^tunnel 



lEi_ 
2U„ 



7 > 1 
7 ~ 1 
7 < 1 



where 

MPI 

TI 

FI 



(3) 



In Eq. (3), Up = 9.33 x lO" 14 IX 2 is the pondermotive 
potential of the laser field, with intensity / in Wcm -2 
and laser wavelength A in /im, and Ei is the ionization 
potential of the atom in eV. The accepted definition is 
that when 7 > 1, the frequency of the laser is greater 
than the tunnelling frequency, hence MPI results. When 
7 ~ 1, the frequency of the laser field is comparable to 
the tunnel frequency of the electron, hence tunnel ion- 
ization is the most prominent mechanism. Conversely, 
when 7 < 1, the laser field comparable to the Coulomb 
field between the nucleus and the electron, therefore clas- 
sical FI dominates. To determine how well the Keldysh 



adiabaticity parameter quantifies ultrafast ionization of 
xenon, our quantification of the ionization mechanism as 
presented in Fig. 4 is converted to adiabaticity in Fig. 5, 
and the data normalized to unity at 7 = 0.1. 

In Fig. Eta), the PPI to Xe™+ (n = 1, 2, 3) in a circu- 
larly polarized laser pulse as a function of 7 is presented. 
Interestingly, the transition between MPI and TI for n — 
1 and n = 2 occurs at the same 7 = 2, thus the identifica- 
tion of the presence of MPI at low intensities in Xe 2+ is 
realistic. The vertical dashed line indicates the transition 
between MPI and TI occurring at 7 = 2, which is consis- 
tent with the intensities at which the transition between 
MPI and TI is observed for Xe+ and Xe 2+ (see Fig. 4). 
The position of this transition is consistent with the pre- 
diction of Keldysh. Furthermore, the PPI of Xe™ + (n = 

1, 2, 3) all tend to unity at 7 = 0.5. For 7 < 0.5, FI dom- 
inates for (n = 1, 2, 3). Fig. 5(a) also contains a 5-point 
running average through the experimental points (solid 
lines) , to indicate the trend of the data. These trends are 
duplicated onto Fig. 5(b) as dashed lines, allowing the 
comparison of PPI in circular and linear polarized laser 
pulses. 

As is apparent from Fig.JSfb), the PPI of Xe™ + (n = 1, 

2, 3) in a linearly polarized laser pulse as a function of 7 is 
rather different from the circular C&1S6, (IS expected from 
Fig. 3 and 4 due to the presence of RI. When n = 1, there 
is little deviation between linear and circular over the full 
range of 7 recorded. However, for n — 2, around 7 = 1.2, 
the PPI breaks away from the response observed in the 
circular case (long dashed line), and tends to follow the 
PPI response for n = 1. This transition, due to RI, has 
not been observed in this manner before. Importantly, 
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the n = 1 PPI is the upper limit of PPI enhancement by 
RI for n = 2. As is apparent from Fig. 5(b), the PPI 
as a function of 7 for n = 3 exhibits a similar behaviour, 
however as there are now two RI processes present, as 7 
is decreased, the PPI for n — 3 tends first to the PPI for 
n = 2 around 7 = 1.3, then around 7 = 2.3, tends to the 
n = 1 response. 



III. SUMMARY 

A new investigation into atomic ionization dynam- 
ics has been presented, employing a novel experimental 
technique allowing a direct comparison between linear 
and circular polarizations. The key to effective inten- 
sity matching (EIM) is defining the spatial dependence 
of ionization yield as independent of polarization type. 
Intensity selective scanning (ISS) is used to measure the 
ionization of xenon as a function of laser intensity by 
translating the focusing optic with respect to a spatially 
limited time-of-flight spectrometer. 

By deriving a solution for the diffraction of a Gaus- 
sian laser pulse through an arbitrary ABCD optical sys- 
tem, geometric effects have been removed from the EIM- 
ISS results, producing partial probabilities of ionization 
(PPI) for Xe n+ (n = 1, 2, 3) for both linear and circular 
polarizations. This technique has allowed clear measure- 
ment of the PPI due to recollision ionization, which con- 
tributes significantly to double and triple PPI in a linear 
polarized laser field. Multiphoton, tunnel and field ion- 
ization contributions are clearly identified for the charge 
states presented for both linear and circular polariza- 
tions. We find that for the lowest intensities, multipho- 
ton ionization successfully predicts ionization to Xe + and 
Xe 2+ . As the laser intensity is increased, tunnel ioniza- 
tion theory is extremely successful, up to the intensity at 
which classical field ionization dominates. 

The volume-independent PPIs have also been pre- 
sented as a function of the Keldysh adiabaticity parame- 
ter for linear and circular polarizations. The applicability 
of multiphoton and tunnel ionization theory allows a pre- 
cise definition of the transition between ionization mech- 
anisms, which occurs when the Keldysh adiabaticity pa- 
rameter, 7 = 2. Furthermore, an interesting dependence 
on the lower ionization states is observed for recollision 
ionization in a linearly polarized laser pulse. This due to 
the necessity to generate the source ion before recollision 
can proceed. 
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APPENDIX A: DIFFRACTION 
QUANTIFICATION 

The conventional approach for treating the propaga- 
tion of light beams is through the solution of the gen- 
eralized Huygens-Fresnel integral (see for example [34|). 
Specifically, the Collins form [35$ of this integral is often 
used, as it allows the ABCD system matrix to be directly 
incorporated into the calculation. When the propagation 
of the beam is limited by some form of aperture, solving 
the Huygens-Fresnel integral analytically becomes more 
difficult, and it often requires a numerical solution. The 
difficulty then becomes one of generating a sufficiently ac- 
curate, computationally efficient solution. This has been 
addressed recently by Lii and Ji [l4| , who deftly adapted 
an approximation of the aperture function of Wen and 
Breazeale 0] to allow the analytical treatment of multi- 
apertured ABCD optical systems for the first time. How- 
ever, the approximate aperture function employed intro- 
duces oscillations into the light amplitude even before any 
diffraction effects are present, thus it is felt that a more 
accurate solution to the problem of propagating Gaus- 
sian beams through apertured ABCD optical systems is 
warranted. 



1. General solution 

In considering the Huygens-Fresnel integral |3jJ for an 
arbitrary ABCD system matrix in polar co-ordinates, the 
beam intensity on the output plane U (r 2 , (p, z) is defined 
as 



U(r 2 ,(p,z) 



1 

AB/o 



U Q {r 1 ,6,z = 0) exp(ifcS') n dn d9 



also referred to as the Collins diffraction integral \3l 
where 



2rir 2 cos[9-<f>]+Dr%) (A2) 



is the path between point (n , 0) on the input plane and 
point (r 2 , <j>) on the output plane, and z is the distance 
between the input and output planes. Variables A, B and 
D are the elements of the system matrix. For a Gaussian 
input beam, in the plane defined by r\ and z — 0, we 
define 



tf o (ri,0,«=O)=exp(-r?/w?) 



(A3) 



where to g defines the beam radius in the input plane. By 
substituting Eq. (A3) and Eq. (A2) into Eq. (Al), and 
employing the following expression from laser resonator 
calculations 36] 



exp 



kr x r 2 
X 



■ cos 



-6) -16 



dB 



2ir\ l exp(— i l<j>) Ji 



krrr£ 
X 



(A4) 
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where J; is a Bessel function of the first kind and I th 
order, we arrive at the following general expression for I 
= 0: 

U(r 2 , </>, z ) = Yb Cxp 1 25 j cxp ' ^ 



/•OO 

/ exp(Pr?) J (Qri)ridri 
Jo 



where 



_J_ ifc^ 
w 2 ^ 2B 



Q 



kr 2 
~~B 



(A5) 



and k = 2tt/\. The analytical solution of this integral is 
now examined for the unapertured and apertured cases. 
It is assumed that the ABCD optical system generates 
a focus in the vicinity of the output plane, however the 
solutions presented apply for any system matrix. 



2. Unapertured solutions 

Many contemporary laser-dilute matter experiments 
assume an unapertured Gaussian pulse profile, which 
when focused generates the following well-known distri- 
bution of intensity, I(r 2 , z) 



In 



i + (z/ zo r 



exp 



-2r 2 2 



W g [h/ [I + (Z/Z ) 



(A6) 



where r 2 and z define the co-ordinate frame in the vicin- 
ity of the focus, and the the beam waist, uj = 2/A/7rw s 
and Rayleigh range, z a — ttljq/X characterize the inten- 
sity distribution. In these expressions, / is focal length 
of the lens generating the focus, A the wavelength and uj g 
the 1/e 2 unfocussed beam diameter. The on axis distri- 
bution can be trivially found by setting r 2 = 0. Solving 
Eq. (A5) for r 2 = 0, the on-axis unapertured solution for 
the current treatment is found by evaluating the integral 

Jo°° dri ' S ivin g 



3. Apertured solutions 

To allow for the action of a finite aperture before the 
ABCD optical system, the integral in Eq. (A5) is evalu- 
ated J °dri where a is the aperture radius. The on-axis 
solution is found by setting r 2 = 0, producing the follow- 
ing expression 



U(r 2 = 0, z) 



XBP 



exp(i/cz)(exp(a 2 P) - 1) (A9) 



where the conditions applied to Eq. (A5) are applicable 
here. 

To find a solution to the off-axis apertured system, the 
integral in Eq. (A5) is evaluated Jgdri where a is the 
aperture radius. This equation cannot be solved analyti- 
cally, so two methods of calculating the spatial intensity 
distribution are employed. Firstly, we perform a Taylor 
Series expansion around r 2 = 0, of the integral term in 
Eq. (A5), so the general form is thus 



U(r 2 ,z) = — cxp 

oo y 



+ Z 



2B 



y=0x=0 VVy ' y ' 



(A10) 



where x and y are the indices of the expansion terms. 
This solution can be used to calculate the distribution 
of intensity for certain optical systems depending on P, 
Q and a. For a relatively sharp focus (/ of the order of 
ujg) with little diffraction (a ^> w 9 ), the number of terms 
required in the Taylor Series is small. However, when the 
system exceeds these limits, particularly when there is 
significant diffraction of the incoming beam, the number 
of terms required becomes untenable, and the (a;!) 2 term 
requires impracticable computational power to evaluate. 
In this case, the intensity distribution may be evaluated 
by converting to a finite element problem through 



U(r 2 = 0, z) 



XBP 



exp(ifcz) 



(A7) 



This expression may be expanded to generate the off- 
axis unapertured solution, by performing the integration 
in Eq. (A5) as before for J^dri, but with r 2 ^ 0, with 
a solution 



U(r 2 ,z) 



2tt 

Jb 



exp 



ifc 



2B 



+ z 



r 1 exp(Pr 2 )J (gr 1 )Ar 1 (All) 



n=0 



U(r 2 ,z) = - 



XBP 



[kDr\ 
2B 



exp j o M " | exp(ifcz) exp(Q 2 /4P) 

(A8) 

This solution requires the real part of P to be negative, 
and the imaginary part of Q to be zero. As this is sat- 
isfied for all cases by the definitions in Eq. (5), this is a 
universal solution. 



This solution requires accurate and rapid evaluation of 
the Bessel Function, which becomes processor-intensive 
as Qri becomes large. By approximating J n {Qr\) ~ 
(2/71-Qri) 1 / 2 cos[Qri - (rwr/2) - (tt/4)] when Qr\ > 8, 
evaluation of this sum becomes very efficient while re- 
taining a deviation from the expected value of < 1 in 
10 4 . The off- axis solution is therefore 
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U{r 2 ,z) = — exp 



2B 



+ z 



„2\ ( 



\7rri( 



1/2 



cos ( Qn - - ) Ari 



(A12) 



r 



By making An small, Eq. (All) and (A12) will produce 
an extremely accurate quantification of the output plane. 
The selection of which form to employ depends on the ge- 
ometry of the ABCD system and the size of the incoming 
beam relative to the aperture diameter. Computational 
efficiency may be further improved through the applica- 
tion of Adaptive Mesh Refinement [33, where the grid 
spacing depends on the rate of variation in U resolved 
along both the r 2 and z axes. 



An important test for the off-axis apertured solution, 
Eq. (A10) to Eq. (A12) is to check for convergence to the 
off-axis unapertured solution, Eq. (A8) and the accepted 
definition, Eq. (A6) when a 3> io g . Through rigorous 
two-dimensional comparisons for a variety of ABCD sys- 
tems, beam radius tu g and aperture radius a have shown 
that all solutions presented are self-consistent, and, when 
applicable, consistent with the accepted solution. 
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